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Effects  of  Boundary  Conditions  and  Initial 
Out-of-Roundness  on  the  Strength  of 
Thin- Walled  Cylinders  Subject  to 
External  Hydrostatic  Pressure 

By  G.  D.  GALLETLY*  and  R,  BART* 


Lsing  classical  *  mill -deflection  theory,  in  investigation 
was  made  of  the  effect#  of  boundary  condition#  and  initial 
out  -cf-rcundne##  on  the  strength  at  cylinder#  subject  to 

in  this  paper  for  initially  cut-cf-round  cylinder#  with 
clamped  end#,  and  a  slifhtly  modified  form  of  the  equa¬ 
tion#  previously  derived  by  Bodner  and  Berk#  lor  iimply 
supported  ends,  were  applied  to  sonic  actual  test  result#  ob¬ 
tained  from  nine  steel  cylinder#  which  hod  been  subjected 
to  external  bydrortatic  pressure.  Three  rcmiempincal 
methods  for  determining  the  initial  out-of-rcundncs*  of 
zhe  vyfinde.'salto  were  investigated  and  these  are  described 
in  the  paper.  The  investigation  indicates  that  if  the 
initial  out-cf-roundness  is  determined  in  a  manner  similar 
to  that  suggested  by  Holt  then  the  correlation  between 
the  experimental  and  theoretical  results  is  quite  good. 
The  investigation  also  indicates  that  while  the  difference 
m  collapse  pressures  for  clsmped -^nd  and  simply  sup- 
jxarted  perfect  w/lindcrs  may  be  quite  considerable,  this 
does  not  appear  to  be  the  case  when  initial  out-of-round- 
ntucs  of  a  practical  magnitude  are  considered. 

Introduction 

SEVERAL  analyses  have  appeared  in  the  literature  for  the 
rlasut  bu-.kJu»g  of  a  thin  cylindrical  shell  subject  vo  external 
hydrostatic  pressure  Cl)  *  The  majority  of  these  analyse# 
have  been  based  upon  the  classical  small-deformation  theory  of 
thm  shells  and  have  assumed  a  geometrically  perfect,  stress-free 
structure  poor  to  loading.  The  correlation  obtained  between 
these  tneories  and  experimental  results  has  been  good  for  long 
cylinders  but  rather  poor  for  thort  cylinders.  Effort#  are  cur¬ 
rently  being  made  by  several  investigators  to  explain  the  dis¬ 
crepancy  between  theory  and  experiment  in  the  short-cylinder 
range  by  the  use  of  large-deflection  theory  At  this  date,  how¬ 
ever,  it  is  not  known  by  bow  much  this  discrepancy  will  be  re- 
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duced  as  nil  the  final  report*  on  their  work  have  not  been  pub* 
lished. 

One  poewble  cause  for  the  di-jcrepancy  botwecu  theoretical  and 

□c#s  of  the  cylinders,  and  a  number  of  investigations,  using  small- 
deflection  theory,  already  have  been  nude  on  the  effect  of  initial 
irregularities  on  the  collapse  pressure  of  cylinders  subject  to  ex¬ 
ternal  hydrostatic  pressure  (2, 3, 4).  As  eras  to  be  expected,  these 
analyses  showed  that  the  initial  irregularities  reduced  the  failure 
pressures  below  those  of  the  perfect  cylinders.  However,  when 
these  analyses  were  applied  to  some  models  which  had  been  tested 
experimentally,  they  predicted  failure  pressures  which  were  less 
than  three  quarters  of  those  observed  experimentally.  Since 
these  analyses  had  assumed  simple  supports  at  the  ends  of  the 
cylinders  and  it  was  probable  that  the  boundary  conditions  of  the 
models  were  somewhere  between  the  extremes  of  simple  supports 
and  clamped  ends,  it  was  of  interest  to  investigate  the  reduction 
in  collapse  pressure  of  damped-end  cylinders  due  to  initial  irregu¬ 
larities,  tc  sco  if  the  assumed  boundary  conditions  rculd  explain 
the  discrepancy  between  experiment  and  theory.  Also,  the 
Analyses  assume  that  the  initial  outrof-roundness  in  the  cylinder# 
is  similar  to  one  of  the  muL*  into  which  a  perfectly  cinmlir 
cylinder  of  the  same  dimensions  would  buckle,  and  actual  shells 
never  satisfy  this  condition.  It  thus  seemed  desirable  to  In¬ 
vestigate  the  various  simplified  methods  that  have  been  sug¬ 
gested  for  determining  the  initial  oubof-roundness  of  the  cylinders 
to  stc  wuat  effect  these  had  upon  the  computed  failure  pressure. 
These  methods  are  described  in  the  paper. 

It  is  also  of  interest  to  note  that  there  are  other  limitations  in 
the  existing  linearised  theories.  These  are: 

(а)  The  fact  that  in  these  problems  there  usually  exist  other 
buckling  pressure#  close  to  the  minimum  buckling  pressure. 
Thus,  use  ot  only  that  mode  of  initial  out-cf-roundnws  which 
corresponds  to  the  minimum  buckling  pressure  is  really  only  de¬ 
fensible  at  pressures  very  ckxe  to  the  minimum  buckling  pressure. 

(б)  The  simple  yield  criterion  used  to  predict  failure. 

This  point  is  discussed  under  the  sectior  entitled  1  Assumptions 
Made  in  Analysis.” 

The  authors  have  not  investigated  the  foregoing  factors  but 
hope  to  do  so  la  the  future.  % 

The  approach  used  in  this  paper  is  similar  to  that  of  Bodner  and 
Berks  (3),  except  that  instead  of  using  a  Donnell-type  equation 
Galerkin’s  method  was  eu  ployed  in  conjunction  with  a  modified 
Donnell-type  eqjition.  The  initial  out-of-roundneas  pattern 
assumed  by  Bodmer  and  Berks  was  of  (he  form 

u*  ••  «  vin  m  $  cos . 

L 

(origin  st  mid-length  j  while  that  assumed  by  the  present  authors 

was 


*  .  flr  2rz"i 
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iORgia  a;  one  end  of  the  cylinder)  Thus,  in  bote  cases,  the 
initial  o'jt-of-roundces?  satisfied  the  relevant  boundary  condi¬ 
tions  and  ni  also  similar  in  form  to  ie  of  the  assumed  buckling 
mode*.  The  magnitude  of  the  initial  out-of-roundnees  at  the  ends 
and  center  of  the  cylinder  was  also  the  same  in  both  cases. 

The  Bcdner»Berks  solution  and  that  presented  herein  thus 
represent  lower  and  upper  bound*  for  the  effect  of  initial  eccen¬ 
tricities  on  the  collapse  pressures  of  elastically  supported  un- 
xtiffeaed  cylinder*,  when  the  initial  eccentricities  have  the  same 
thape  as  one  of  the  assumed  buckling  inodes  of  the  perfect 
cylinder.  'While  the  two  solutions  are  not  exact,  they  should 
provide  good  approximation*  to  the  exact  solutions.  One  of  the 
limitation*  of  using  Donnell's  equation  is  that  the  number  of  cir¬ 
cumferential  lobe*  should  be  fairly  high,  and  thus  the  results  will 
be  shghlly  in  error  for  very  long  cylinders  which  buckle  into  two 
or  three  circumferential  lobes. 

The  final  results  of  the  investigation  axe  given  m  Fig  2  and  in 
Tables  2,  4,  and  5.  It  can  be  seen  that  the  correlation  between 
experiment  and  theory'  is  quite  good  when  method  (c)  is  used  to 
determine  the  initial  out-of-roundness  of  the  models.  (Method 
(e)  is  similar  to  that  suggested  by  Holt  (5).J  However,  it  is  not 
claimed  that  the  results  give  a  complete  answer  to  the  problem 
and  more  work  of  both  an  experimental  and  theoretical  nature  is 
required 

Method  or  Analysis 

The  modifications  to  Donnell’s  equation  brought  ebowt  by 
initial  eccentricities  in  the  shell  have  been  presented  by  Bodner 
and  Berks  (3X  and.  prior  to  them,  by  Cicala  (6).  The  equations 
also  have  been  derived  by  the  authors  in  the  Appendix,  using  a 
somewhat  different  approach  to  that  adopted  by  the  previously 
mentioned  author*  Tor  the  case  of  uniform  external  hydrostatic 
pressure  applied  on  ail  *id«  of  an  imperfect  cylinder  the  relevant 
equations  are,  from  Equations  (19c),  (196},  (21],  and  (24)  in  the 
Appendix 
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h  *•  thickness  of  ebefi 
R  -  mean  radius  of  shell 
p  «  applied  hydrostatic  pressure 
E  -  modulus  of  elasticity 
X>  -  Eh*/ 12(1 —»») 
v  •»  Poisson's  ratio 

F  -  stress  function  of  the  total  membrane  stresses 
/  s»  i  a«\* 

\j>.  m  initial  radial  out*of-rouodt»ew  (+  inwards) 


u,  v,  v  «  elastic  axial,  tangential,  and  radial  (+  inwards) 
displacements  of  the  imperfect  cylin ,  :  minua  the 
uniform  compression  experienced  by  »  perfect 
cylinder  (eso  Equation  (IS)  in  the  Appendix). 

The  t  scripts  x  and  0  indicate  partial  differentiations  with  re¬ 
spect  to  those  variables. 

The  patterns  assumed  for  w  and  u*,  and  which  satisfy  the 
boundary  conditions  for  clamped -end  cylinder*,  were  as  follows 


w*  — -  rin  md 


B  -  half  amplitude  of  w-displacement 
e  -  maximum  value  of  initial  radial  out-of-roundness 
m  «■  number  of  circumferential  waves 
L  »  Unsupported  length  of  shell 
x,  6  -  axial  and  angular  co-ordinate# 

If  the  Expressions  (2)  happen  to  be  an  exact  solution  of  the  prob¬ 
lem,  then  they  will  satisfy  the  differential  equation  of  equilibrium, 
Equation  lie),  exactiy.  However,  as  both  tc  and  ♦?*  were  chosen 
to  satisfy  the  boundary  conditions  rather  than  the  equilibrium 
equatioo,  this,  in  general,  will  not  be  the  case.  The  resulting  ex¬ 
pression  will  be  a  function  of  x  and  C  which  we  shall  denote  by  Q. 
Galerkin's  equation  for  determining  the  relations  between  the  co¬ 
efficients  B  and  t  is  then 

f1'  fLQ* in  ^[i  —  C0«  — j  RdMx  -  0 . 13) 

where  i  assumes  the  value*  1, 2, 3 . 

for  i  *  m  Equation  (31  will  be  found  to  be  zero  identically 
For  i  -  m  the  following  relation  between  B  and  e,  obtained  from 
liquation  (3J,  will  Ire  found  to  hold 
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where  p.t «  given  by  the  expreeeion 
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The  smallest  value  ol  the  buckling  pressure  of  the  perfect 
cvlinder  is  found  by  minimizing  ’“•\uation  (6)  with  respect  to 
m.  A  relation  steiibr  to  that  expressed  by  Equation  (6)  ha*  re¬ 
cently  bf«i  presented  by  Nash  (7)  using  an  energy  method.  A  re¬ 
lation  similar  to  Equation  |4 )  vas  also  obtained  by  Bodner  and 
rlcrks  f::  rjr.ply  suppcrb*d  imjx'rfect  cylinder*. 

Tlius,  from  Equation;  {2)  and  [4],  we  obtain  the  following  ex¬ 
pression  for  ta 

v>  ••  ~  —  — —  rin  m  6  f  1  —  cos  ~  I.  (5) 
i  i>*  —  P  L  h  J 

The  bending  moments  in  the  shell  c as  then  be  calculated  from 
the  relations 
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The  muimnm  bending  stresses  are  then  given  by 
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lo  obtain  the  total  norail  ftresaee  t re  now  uave  to  add  the 
membrane  nrew*  to  Equation  {8j.  To  devenmne  these  Utter  we 
solve  Equations  [idj  and  [Sj  for  the  stress  funcUon  P  (periodic 
tens*  only;*  The  total  membrane  stresses  arc  then  given  by 


<r.. 
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The  total  normal  itreaoee  are  obtained  by  adding  algebraicall> 
Equations  (Sj  and  (9j.  The  greatest  normal  streeaes  occur  at 
cud-length  of  the  cylinder  (x  -  1/2)  and  where  tin  md  •  ±1 
v trough  and  crest  point*  of  the  lobe*>.  At  these  point*  the  twist¬ 
ing  moment  is  rero  and  thus  the  normal  stresses  are  principal 
stresses.  The  absolute  maximum  normal  stresses  occur  at  the 
outer  shell  wall  for  the  trough  point*. 

Having  obtained  the  maximum  principal  etresae*  o,  and  <r#  m 
terms  of  the  initial  out-of-roundnees,  tie  geome-nc  parameters 
and  the  applied  external  hydrostatic  p,’  essure,  we  now  employ  the 
octahedral  shear-stress  criterion  of  failure  (which  give*  the  same 
result*  a*  the  Hencky-vcn  Mitel  criterion  of  failure),  vis. 


<rt*  •*  <r#*  +  s,* —  vt&i- •  . (101 

where  <rf  1*  the  yield  pemt  of  the  material.  Substitution  of  the 
ffitnaian  principal  stresses  om  and  0%  in  Equation  (10 j  then  gives 
an  equation  relating  the  initial  out-of-roundntts,  the  geometric 
parameter*  of  the  shell,  the  yield  point  of  the  material  and  the 
pressure  at  which  the  she'd  begins  to  yield  pw. 


It  should  be  noted  that  instead  of  uring  the  yielding  criterion 
given  by  Equation  (10)  where  am  and  c>  are  principal  afxessea,  it 
U  more  accurate  to  use  the  expression 

<V  *  +  s'.*  —  <rff(r,  v  3r^» . [IOoJ 

where  now  <r„  a#,  and  r,#  are  the  normal  and  shear  stresses  at  any 
point  and  which  are  functions  of  z  and  6.  Yielding  will  dirt  occur 
in  the  cylinder  for  those  value*  of  z  and  0  which  maximize  the 
right-hand  eide  of  Equation  [10a].  However,  to  compute  these 
value*  of  x  and  6  by  differentiation  of  equation  {10a  |  involves 
more  complications  than  seem  warranted.  Trials  indicate  that  the 
stress  condition  at  the  outer  shell  wall  for  tiough  points  of  a  lobe  Is 
probably  aa  unfavorable  a*  anywhere  else.  As  mentioned  earlier 
the  twisting  moment  is  zero  at  these  points  and  thus  Equa¬ 
tion  [10a)  reduces  to  Equation  (10). 

As  we  shall  later  present  curves  of  pr,  the  pressure  at  which  shell 
yielding  commences,  versus  */h,  the  initial  eocentridty- 
shetl  thickness  ratio,  for  both  simply  supported  and  damped-end 
cylinders,  we  have  summarised  the  results  obtained  in  this  paper 
and  those  obtained  by  Bodner  and  Berks  in  Tablo  1.  (We  have 
added  e  few  terms  to  the  latter  solution,  aa  Bodner:  and  Berks 
neglected  the  periodic  terms  in  Equation  [9]).  Alio,  we  used 
the  expression  t c^fR*  for  the  circumferential  change  in  curva¬ 
ture  instead  of  («>*/#*  +  r t/R)  which  was  used  by  Bodner 
and  Berks  when  computing  the  bending  stresses  due  to  initial 
out-of-roundness.  Our  expression  is  consistent  with  the  expres¬ 
sion  used  in  deriving  the  approximate  equations  of  equilibrium 
for  an  initiaily  out-of-round  cylinder  (eee  Appendix)  aod  m«o  is  the 
same  as  that  used  by  Donnell  in  his  work  on  perfect  cylinders  (8). 
The  effect  of  using  vh/R*  instead  of  (wh/R*  +  tt/R)  is  to 
eliminate  the  quantity  f  which  appears  in  the  equations  de¬ 
veloped  by  Bodner  and  Berks.)  It  can  be  seen  that  the  final 
equation  (Equation  (12!  in  Table  1)  relating  the  initial  out-of- 
roundness  and  the  pressure  to  cruse  first  yielding  is  essentially 
the  same  for  both  clamped-end  and  simply  supported  cylinders, 
except  for  the  factor  of  2  required  by  the  definition  of  «  as  the 
maximum  initial  vut-of-roundness  and  tho  slight  changes  in 
definition  of  the  quantities  K  and  0  appearing  in  that  equation. 
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Besides  tbs  Approximate  strain  and  change  in  cur**ture-dia- 
plissmsnt  relations  used  in  the  derivation  of  the  Donnell  equa¬ 
tion  for  an  initially  out-of-round  cylinder,  and  the  assumption 
that  it  ij  initially  street  free,  it  might  be  useful  to  lint  some  of  the  * 
othci  assumptions  that  have  been  made  both  in  the  analysis  of 
this  paper  and  that  of  Bodnar  and  Berks.  Thsoe  are  as  follows; 

1  That  the  circumferential  membrane  stress  in  the  shell  is 
constant  along  ita  length  and  equsl  to  — pR/A,  whereas  this  is 
actually  not  the  ease. 

2  The  assumption  that  the  initial  out-of-roundnees  is  small 
(with  an  order  of  magnitude  of  one  shell  thickness)  if  symmetric 
with  respect  to  the  center  line  of  the  shell,  and  has  the  same  form 
ae  one  of  the  buckling  modco  of  a  perfect  cylinder  with  the  tame 
shell  dimenrions.  Actual  rhslls  rarely,  if  ever,  satisfy  the  last  two 
requirement*  and  so  the  question  arises  as  to  how  the  initial  out- 
of-roundnes*  should  be  waasored 

3  The  assumption  that  failure  occurs  (appearance  of  risible 
lobes)  when  the  most  highly  stressed  points  in  the  cylinder  start 
to  yield.  Actually,  failure  does  not  occur  until  plastic  regions 
form  at  the  trough  and  crest  points  of  the  lobes.  The  preesure  re¬ 
quired  to  produce  these  yield  soces  is  greater  than  that  at  which 
the  most  highly  stressed  points  begin  to  yield  and  neglect  of  this 
effect  therefore  underestimates  the  strength  of  the  shells.  An 
adequate  theory  to  take  this  effect  into  account  has  not  been  de¬ 
veloped  as  yet,  but,  as  for  beams,  presumably  the  ratio  of  the 
preanozv  to  cause  first  yielding  to  the  p.cssure  required  for  the 
formation  of  plastic  regions  depends  on  the  relative  magnitudes 
d  the  direct  stresses  m  the  cylinder  wall  and  the  bending  stresses 
resulting  from  initial  out-of-roundnea*. 

4  The  assumption  that  Poisson's  ratio  is  a  constant  and  equals 
0.3. 

As  zt  is  intended  to  apply  the  analyses  developed  for  unstiffened 
cylinders  to  stiffened  cylinder*  which  failed  by  buckling  between 
nng  stiffeners,  it  would  be  well  if  we  enumerated  the  additional 
assumptions  that  were  made.  These  are: 

5  That  the  stiffening  rings  at  the  ends  of  any  bay  are  per¬ 
fectly  circular,  ae.,  they  do  not  have  any  initial  out-of-roundnes*. 
This  never  occurs  w  practice,  of  course,  but  should  not  be  too 
eenous  if  the  circularity  of  the  stiffening  rings  u  very  much  better 
than  that  of  the  shell,  or  if  'he  predominant  mode  of  initial  o"t-of- 
roundness  in  the  nngs  la  very  different  from  the  predominant 
mode  in  the  shell. 

5  As  for  unstiff cned  cylinders,  the  drcumfertnUsl  membrane 
strtss  in  the  perfect  cylinder  is  assumed  to  be  — pR/h,  wnereaa 
it  actually  vanes  along  the  length  of  the  shell.  A  more  correct 
repraentrnua  of  the  streea  distribution  wouid  be  obtained  b> 
ucy  ihe  analysis  of  too  Banden  and  CKlnther  (0>,  or  more  ac¬ 
curately  ftiU,  that  of  Salerno  and  Pulos  (10). 

Mztbods  or  DrrraKiM.NO  Ixitlj,  OuT-or-RotNOsaas 

Az  mentioned  hitherto,  the  analyses  assume  that  the  initial 
out-05 -roundn^e*  is  symmetrical  about  the  mid-length  of  the  shell 
:zji  that  itc  circumferential  variation  b  in  the  shape  of  one  of  the 
buck’ing  modes  of  a  perfect  cylinder.  Actual  sheds  do  not  meet 
either  of  the**  requirements  If  we  do  not  make  a  harm  otic 
aatlyxs  J  Ox  initial  out-of-rpundneas,  and  abo  extend  the  theory 
to  account  for  the  various  harmonic  components,  the  qui^Uou 
ariaes  as  \0  how  we  shall  measure  the  quantity  e,  which  la  defined 
in  the  analysisas  the  maximum  initial  emt-of-roundnes*  when  its 
shape  is  to  ons  of  the  buckling  modes.  As  far  as  the 

authors  are  aware  three  simplified,  semicmpirlcal  methods  foi  de¬ 
termining  the  initial  out-of-roundnea*  have  been  p.v posed  »«  the 
literature  so  far.  These  are  (ace  Fig.  1): 

The  centroid  of  the  rnituu  circuomy  <x-*toui  is  first  deter- 
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METHOD  (0) 

Fio  1  iLLCtnuTxoxs  or  Tbsxx  Mrracns  roa  Ditxjucnino 
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BR-6 

mined  Then  the  angle  r/m,  where  m  is  the  number  of  iobes  into 
which  the  perfect  cylinder  would  buckle,  if  calculated.  A  sector 
of  a  circle,  subtending  this  angle  r/m  is  then  drawn  on  transparent 
paper  and  placed  with  its  apex  at  the  centroid  of  the  initial  circu¬ 
larity  contour.  The  sector  is  then  rotated  so  that  it  traverses  the 
entire  circumference  of  the  circularity  contour  until  the  location  Is 
found  at  which  the  maximum  difference  between  the  two  sector 
radii  occurs.  The  initial  out-of-roundnea#  is  then  taken  as  this 
maximum  difference.  This  method  of  determining  the  initial 
cut'cf-roundneas  is  essentially  that  proposed  by  Saunders,  Tril¬ 
ling,  and  Windenburg  (11, 12). 

(6)  Both  the  centroid  and  the  area  of  the  initial  circularity  con¬ 
tour  are  determined.  The  radius  Rm  of  the  circle  whoso 
area  is  the  same  as  that  of  the  initial  circularity  contour 
is  then  determined.  A  circle,  with  center  at  the  centroid  of 
the  initial  circularity  contour  and  of  radius  ft*,  is  then  drawn. 
The  initial  out-of-roundneaa  is  then  taken  as  the  maximum  value 
of  (/?,  —  /?*)  and  (ft*  —  A,  J,  where  R,  and  Rt  are  the  radius  vec¬ 
tors  from  the  centroid  to  points  on  the  initial  circularity  contour 
which  are  exterior  and  Interior,  respectively,  to  the  circle  of  radius 
Rm  A  method  for  determining  the  initial  out-of-rcundneas  simi¬ 
lar  to  the  foregoing  has  been  suggested,  among  others,  by  Bodner 
and  Berks  (3). 

(e)  As  in  (6)  both  the  centroid  of  the  initial  circularity  contour 
and  the  radius  Ru  of  the  mean  circle  are  determined.  Abstain 
(a),  the  angle  r/m  is  calculated.  The  arc  length  of  one  halMobe 
is  then  obtained  aa  ( r/m)Rm ,  This  arc  is  then  moved  around  tho 
initial  circularity  contour  with  its  end  points  always  In  contact 
with  the  contour.  The  initial  out-of-roundneas  la  then  taker,  as 
the  maximum  radial  distance  between  the  circularity  contour  and 
the  arc.  This  method  for  determining  the  initial  out-of -round ncas 
is  somewhat  similar  to  the  method  proposed  by  Holt  (5). 

Numerical  Rxaui/ns 

In  this  section  we  present  the  results  rbtubud  by  applying  the 
anal/tb  developed  for  simply  supported  imperfect  cylinders  (3>, 
and  its  extension  to  clamped  ends,  to  »tsel  wUded  cylinders 
that  have  been  tested  at  the  Taylor  Model  Butin.  At  failure,  all 
these  models  had  lobes  which  partially  covered  the  dreumferenoe 
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»  No  datiuction  suuie  m  Wiadaabur*  •  aotM  between  buckUai  end  audaua  ruiUlntag  prtcrura  for  theoe  taodtij,  fceaoe,  they  were  usuxatd  (9  b*  eqcii 


Cfrtculeied  fox  median  rarfeoe  of  eheU  it  mid-bey  ieAjth  aiu4  eoefyei*  of  roa  Send  an  end  Odnthix  la  ©oajuaetton  with  ootehedral  lieu  street  eriwrlot 


TaB»-E  I  MAXIMUM  «,  LVALUES  OBTAINED  BY  DIFFERENT  METH3D8  FOR  DETERMINING  INITIAL  OUT-OF-ROUNDNES8 

(1) 


Model  no. 
SO 
** 

60 

71 

42 

61 

BR-i  re. 

BR-4  [4. 

BR-S  4. 


S)» 

10) 

6.  1C) 


(*> 


(S) 


(4> 


.  «)  . 


<«) 


- - Method 

C.E 

<a)8.8 r 

Method  (5) 

'c.z.® 

Outward  — > 

8-B. 

nri-  ■-» 
8.8. 

• - Mai.  of 

O.E. 

0  1*6* 

0.140 

0.175 

0.016 

0.017 

0.074 

0.07* 

0.074 

0.078 

0.C97 

0.104 

0.09* 

0.057 

0.078 

0.0*6 

0.003 

0.089 

0.093 

0.04* 

0.047 

0  041 

0.041 

0.046 

o.oi> 

0.0*7 

0.041 

0.045 

0  1*9 

0.141 

0.11* 

0.006 

0.011 

0.090 
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0.005 

C.267 

0.272 

0  101 

0.080 

0.001 

0,114 

0.121 

0.114 

0.121 

0  126 

0  1*0 

0.092 

0  027 

0.C2* 

0.C22 

0.027 

0.027 

0.03* 

0  575  (6>* 

0  616(8) 

0  *54(2) 

0.210(6; 

0  *76(6) 

0  267(0) 

0.902(6) 

0.110(6) 

0.07*h) 

0.2*6(10) 

0.376(6) 

0  164(4) 

0  176(4) 

0  765(4) 

0.128(2) 

0.602(4) 

0  069(4, 

0.069(4) 

0.073(4) 

0.210(4) 

0.082(4) 

0.082(4) 

0  741(4} 

0  236  (1C 

)  0  283(10) 

0.247(4) 

0.285(10) 

•  Tabulated  ntu«a  aeecreu  to  epproxiratveiy  ±5  per  oeat.  .. 

i  Nuabera  in  *,  <+r*  bracket*  indicate  the  eutione  it  which  iobee  first  appeared  ia  the  oultibay  i  tiffined  cylinders.  Numbers  in  puentheece  t re  the  ela¬ 
tion*  at  which  the  Eu.xi-'-ia  «/A-rilut*  occurred,  according  to  the  method  ueed. 
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TABLE  4  COMPARISON  OF  PREDICTED  PRESSURES  PSD  FOR  OCCURRENCE  OF  A  VISIBLE  LOBE— USING  THE  DIFFERENT 
- - - - 0UT-0F-R0UNDNE8S— WITH  THE  EXPERIMENTAL  PRESSURES  \P8I) 

Modal 
no. 

SO 

» 

66 
•n 
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G1 
BR-1 
BR-4 
BR-S 
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65 
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68 
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66 

39 

59 
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80 
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57 

50 

69 

61 

92  5 

S 

76 

80 
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o.«o 

0  763 

1.16 
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of  the  shell  The  geometric  ratio*,  yield  points,  experimental  and 
theoretic*:  collapse  pressures  for  these  cylinders  are  given  in 
Table  2  The  first  six  models  in  this  table  were  tested  some  °0 
years  ago  by  Wind en burg  and  Trilling  (11/,  although  they  did 
not  investigate  theoretically  the  effect  of  initial  out-of-roundnes* 
on  the  collapse  pressure  The  last  three  models,  which  are  multi¬ 
bay  cylinders,  have  been  tested  recently  at  the  Taylor  Model 
Basin  (13,  M). 

A  comparison  of  columns  6  and  8  in  Table  2  shows  a  consider¬ 
able  discrepancy  between  the  experimental  and  theoretical 
buckling  pressures  for  even  the  simply  supported  cylinders  For 
some  models,  it  would  also  appear  that  axisymmetnc  yielding 
rather  than  buckling  was  the  cootrolling  mode  of  failure.  This 
can  be  teen  by  comparing  columns  8  and  10  in  Table  2  However, 
it  mil  be  seen  later  when  out-of-roundae«  Is  taken  into  account 
that  the  theoretical  pressures  for  buckling-Jype  failure*  are  lower 
than  the  axisymmetrie  yield  pressures.  It  U  also  of  interest  to 
uote  that  Models  42  and  61  are  the  only  models  for  which  the  ex- 
pcnmenUi  budding  pressures  are  higher  than  those  predicted 
theoretically  for  simply  supported  cylinders,  although  the  sai  > 
supports  were  used  for  these  two  models  as  for  she  four  models 
preceding  them. 

Using  the  geometric  ratios  and  yield  points  ehown  in  Table  2 
and  Equations  [11]  and  |12,  in  Table  1,  it  is  possible  to  construct 
curves  showing  the  relation  between  the  pressure  at  which  yielding 
first  occurs  m  the  cylinder  wall  pt,  and  the  ratio  of  the  initial  out- 
of-roundness  to  the  shell  thickness  «/A.  Two  such  curve*  are 


ehown  in  Fig.  2  for  illustrative  purposes.  In  constructing  these 
curves,  the  value  of  m  ueed  in  Equation  (12]  waa  the  value  of 
m  which  minizrdxed  Equation  [11  ]  in  Table  1.  Theec  values  of  m 
are  Hated  In  parentheses  in  columns  8  and  0  of  Tab'e  2.  These 
values  of  m  do  not  actually  give  the  minimum  prh  a  given  s/A. 
This  point  will  be  dlscuseei  later.  It  also  should  j*  noted  that 
for  t,  h  «•  0  aome  of  the  cum*  in  Fig.  2  do  not  attain  the  elastio 
buckling  pressures  tabulated  in  columns  8  and  9  of  Table  2. 
When  this  is  found  to  occur  it  means  that  the  pressure  to  cause 
axisymmetrio  yielding  of  the  shell  ii  lower  than  the  elastio 
buckling  pressure. 

In  Table  3  are  tabulated  the  maximum  s/A- values  obtained 
using  the  different  methods  for  determining  the  initial  out-rf- 
roundness  described  earlier.  The  «,A~values  were  determined  at 
mid-length  for  the  uastiff eneA  cylinders  and  at  mid-length  of  the 
bays  for  the  multibay  cylinders.  Two  value*  are  listed  fui  each 
model  under  methods  (a)  and  (:)  because,  in  these  methods,  toe 
number  of  lobes  into  which  the  perfect  cylinder  would  buckle  is 
used  and  — ils  number  is  usually  different  for  simply  supported 
and  damped  ends. 

How,  selecting  the  experimental  b  idding  pressures  listed  in 
column  6  of  Table  2  in  conjunction  vnth  their  corresponding  eo- 
centncitie*  listed  in  columns  2, 3,  and  6  of  Table  3,  one  can  plot 
points  on  curves  similar  to  Ffg.  2  which  represent  values  deter¬ 
mined  experimentally .  In  Ta^le  4  we  alto  give  «  numerical  com¬ 
parison  of  the  theoretic*!  and  experimental  pressures  at  which 
risibb  lobes  first  occur.  The  theoretical  value*  in  this  table 


EXPERIMENTALLY  determined  */h  VALUES 


Fio  2  Sjxxll  Yituwo  P»EMuax  toe  Models  33  and  60  Versus 
Eccentricitt*  Ssxix-TaicrxrM  Ratio 

TABLE  5  AVERAGE  VALVES  OF  COLUMNS  («.  (7).  AND  (*)  IN 
TABLE  4.  CLASSIFIED  ACCORDING  TO  CYLINDER  TYPE 
Cyioder  —Method  (o>— -  —Method  (6)—  —Method  <<>--» 

type  C.E.  88  C.E.  8.8  CE  8.8. 

U&aUSrsed  1C*  0  85  l  12  0*9  1  !8  0  01 

Suffered  085  0  76  u  88  079  108  0.91 


were  obtained  from  the  e/h-\  alues  lilted  m  Table  3  in  conjunction 
with  the  theoretical  pw  versus  e/h  curves  similar  to  Fig.  2.  The 
last  three  column*  in  Table  4  snow  the  retioe  of  theoretically  pre¬ 
dicted  pressure*  for  the  occurrence  of  a  viiible  lobe  to  those  ob¬ 
tained  experimentally,  according  to  the  various  method*  used  for 
determining  the  initial  ouLof-rouadness.  The  average  value*  of 
these  last  three  column*,  classified  according  to  whether  the 
cylinder*  were  stiffened  or  not,  are  tabulated  in  Table  5.  It  can 
be  seen  from  Tables  4  and  5  that  use  of  method  (a),  with  the 
assumption  of  simply  supported  ends,  was  the  most  conservative 
in  most  caeca  and  predicted  pressures  which  were  always  below 
thooe  obtained  experimentally.  However,  the  best  correlation 
between  the  experimental  results  and  the  simplified  theories 
discussed  m  this  paper  appears  to  be  obtained  when  method  (<  t  is 
used  for  determining  the  mitral  out-of-roundneas  and  the  <.j  finders 
are  assumed  to  be  simply  supported 
It  was  mentioned  earlier  that  for  a  given  t,  h  the  value  of  m  that 
would  give  the  lowest  pt  was  not  necessarily  the  value  of  m  which 
minimized  the  expression  for  pw  v Equations  tl  1  j  *n  Table  1  >  It 
was  also  noted  at  what  tune,  however,  that  the  error  obtained  by 
assuming  uns  to  he  so  was  nut  vtiy  great  and  aiso  greatly  re¬ 
duced  the  computational  labor.  Some  idea  of  the  *.«or  .nvolved 
can  be  obtained  by  referring  to  Fig  3.  The  curves  m  this  figure 
are  plots  of  m  which  minimize  r,  h  for  a  given  p.  For  a  pt  80  psi 
it  can  be  seen  that  the  minimum  value  of  t,  K  is  0  34  at  m  ••  1$, 
while  at  m  -  16  Ktbe  value  which  minimizes  p«,  the  value  of 
r/Au0  4.  Howevei,  if  we  now  fix  the  value  of  e, h  at  O.t,  then 
tne  mim mu n.  value  of  pr  is  78  ps*  and  occurs  at  m  *  i!#  Similar 


Fio.  3  Initial  EocxxT*icmr-Sro,x.  TmcxNrsa  Ratio  Vex* vs 
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Fio.  4  Ttsjcal  Initial  Circulasitt  Contour*  roa  Model 
DIM 

i  emits  are  obtained  with  the  other  curves  and  also  with  the  m 
versus  e,  h  curves  for  damped-end  cylmders  which  ws  have  not 
included  here.  It  is  thus  seen  that  assuming  the  value  of  m  which 
minimizes  pw  will  also  minimize  pw  for  a  given  t/h  if  slightly  on  the 
unsafe  side,  the  magnitude  of  the  error  depending  on  the  value  of 
t,  h  taken  However,  owing  to  the  fact  that  actual  cylinders  under 
uydrcstaUc  pressure  collapse  in  substantially  the  tame  number  of 
tubes  as  is  predicted  by  theory  for  the  perfect  cylinder,  and  also  for 
simplicity  in  cal  culm  Luna,  we  have  ignored  this  discrepancy. 

It  aUo  will  be  remembered  that  in  applying  the  analyses  to 
stiffened  cylmders  we  made  the  assumption  that  the  stiffening 
rings  were  perfectly  circular.  This,  of  course,  never  occurs  in 
practice  Some  idea  of  th*  degree  of  circularity  actually  present 
can  be  obtained  by  reference  to  F*g .  4.  This  figure  shorn  the 
imual  circularity  oon tours  of  the  stiffening  rinp  bordering,  and  the 


taeil  u  the  center  ot  one  of  t he  bay*  in  which  iobea  firtt  appeared 
m  Model  BR-4  If  we  adopt  method  (e)  described  earlier  u  our 
criterion  for  initial  out-of-roundne».  then  the  stiffening  ring*  of 
this  model  haa  about  one  tenth  the  ou.t-of-rom>dne*s  of  the  shell. 
Similar  results  also  were  obtained  for  the  other  models.  Thus,  for 
the  problem  of  interring  eollapee  of  stiffened  cylmders,  the  as¬ 
sumption  of  zero  initial  out-of-roundne**  o:  the  stiffening  rings 
appears  to  be  a  reasonable  one. 

It  might  be  thought  that  a  haimomc  analysis  of  the  initial  out- 
of-r^undne*5  would  show  that  the  amplitude  of  the  harmonic 
component  corresponding  to  the  value  of  rn  which  minimized  ;v 
(Equation  {1 1 )  in  Table  1 )  was  many  times  that  of  the  other  com¬ 
ponents.  To  investigate  this  point  a  harmonic  analysis,  using  72 
subdivisions  of  the  circumference,  was  made  of  the  initial  circu¬ 
larity  contour  at  Station  4  of  Model  BR-4  uzing  Filon’6  method 
\l 5/  rhe  number  of  waves  into  which  a  perfect  cylinder  wits 
shell  dimf-iSion*  similar  to  BR-4  would  bucsile  is,  accord  mg  to  the 
Imear  theory  used  herein,  11  for  simply  supported  ends  and  12  for 
clamped  ends  From  the  harmonic  an  alms  it  was  found  that 
there  were  nine  narmoxncs  with  amplitudes  greater  than  the 
eleventh  and  two  greater  than  the  twelftn.  It  also  was  found  that 
the  amplitude  of  the  largest  harmonic  (m  -  8)  was  more  than 
twice  that  of  th*  twelfth  mode  and  mere  than  five  times  that  of 
the  eleventh  mode  Furti— rf  even  this  largest  amplitude  was 
much  too  email  to  effect  a  reasonable  correlation  between  theory 
and  experiment  Thus  the  harmonic  analysis  of  the  BR-t  initial 

did  not  produce  any  results  which 
might  be  useful  in  practice. 

In  conclusion,  it  might  be  of  interest  to  mention  that  some  at¬ 
tempts  have  been  made  to  determine  experimental!)  the  longi¬ 
tudinal  form  of  the  buckling  displaeemnte  in  the  multibay  cylin¬ 
ders  The  results  of  these  few  investigations  are  summarized  in 
reference  ( 1 4  However,  more  experimental  wort  still  remains  to 
be  done  before  any  conclusion  reganLng  the  shape  of  the  buckling 
displacement  can  1*  made 

Summary 

Ir.  the  preceding  actions  an  attempt  has  been  made  to  explain 
some  of  the  discrepancies  tba*  exist  between  experimental  and 
theoretical  results  for  cylinders  subjected  to  external  hydrostatic 
pressure  To  do  this  it  was  assumed  that  (c;  the  actual  boundary 
conditions  were  somewhere  between  the  extremes  of  simple  sup¬ 
ports  and  clamped  ends,  (6,  the  initial  out-of-roundness  was 
similar  m  form  to  one  of  the  modes  into  *t»'h  a  perfect  cylinder 
would  buckle,  (cl  the  strew  dit-»ouuon  in  the  equilibrium  prob¬ 
lem  lor  the  perfect  cy  linear  could  be  represented  by  the  membrane 
stresses,  and  (di  tb*  linear  imaiiKiJ'wmatior,  equations  of  equi¬ 
librium  would  describe  the ;  rotlem  adequately  Any  cold-work¬ 
ing,  residual.  or  weldug  stresses,  or  any  elastic  nonhomogeceity 
that  might  have  teen  present  were  neglected  It  also  was  as¬ 
sumed  that  fauure  i formation  of  a  lobe,  would  occur  when  the 
stresses  at  the  most  highly  stressed  point  satisfied  the  octahedral 
shear -strew  entenon  Three  simplified  methods  of  measuring  the 
initial  otst-of-roundnew  were  also  investigated  The  simplified 
analyses,  together  with  the  different  methods  of  measuring  vut-cf- 
roundness,  were  applied  to  mre  steel  welded  cy  finders  with  length- 
diameter  ratios  cf  1  »  to  2,  thrchnesa-diarn'to  ratios  of  0  0025  to 
0  and  '*iei<3  points  of  the  fteei  30,000 1  j  60,000  psi  The 
correlation  between  experimental  and  theoretical  results  was 
quite  good,  w»en  method  .c  w'  used  lor  determining  the  initial 
out-o' -round ness  c*  the  e\fin'  en 
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Appendix 

Affroximate  Equations  or  Equilibrium  for  an  Initially 
Out-of-Round  Cylinder 

These  equations  have  been  derived  previously  by  Bodner  and 
Berks  (3;  ueing  a  different  approach  Our  reason  for  this  new 
derivation  it  to  make  clear  what  approximations  and  assumptions 
are  involved  in  the  simplified  equations  of  equilibrium. 

Consider  a  cylinder  of  mean  radius  ft  which  has  a  small  radial 
iniUal  out-of-roundness  w,  Selecting  x,  0,  and  t  as  co-ordinate 
axe*  and  denoting  the  longitudinal,  tangential,  and  radial  elastic 
displacement*  by  V,  V,  and  W,  it  is  then  possible  to  show  that 
the  strains  at  the  middle  surface  are  given  by 
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I.  -  V.  -r  -J  [17  f  IT.’!  1-  tr,tr>. 

[r  -t-  w-.l  +  ^v  ir,-trj  [  I!3| 
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_  Ll1’*  x  t'->,  A  tr.tft, 

R‘~  ^  R  S'  B  +  R* 

We  shall  amplify  *b«e  complicated  expressions  to  the  following 

«.  -  l ,  r  -j  BV  +  W.tc-, 

«  -  4  +  -^T  +  HJT  1  l>4! 


where  tf,  T,  tad  IP  are  the  displacements  which  would  occur  in 
the  equilibrium  problem  of  a  perfectly  circular  cylinder  under 
uniform  external  pressure 

The  total  potential  of  the  syrtera  UT  i®  then  obtained  by  add¬ 
ing  the  extension&l  and  bendings  train  energies  of  the  shell  and 
subtracting  the  work  done  by  the  external  pressure.  In  calculat¬ 
ing  the  extensional  energy,  we  retained  the  terms  in  u,  r,  to,  tc* 
through  the  second  order  and  the  terms  in  0,  P,  W  (directly  pro¬ 
portional  to  the  applied  pressure)  through  the  firtt  order  only; 
also,  w  j  neglected  the  effect  of  the  deflection  of  the  shell  between 
supports  on  the  displacements  and  stresses  of  the  perfect  cylinder. 

Variation  of  Ur  with  respect  to  u,  v,  and  tc  then  give*  the  dif¬ 
ferential  equations  of  thi  problem.  Further  manipulation  of  these 
equations  results  in  the  follow  ng 

**“ -¥*"•-*?  (a) 

VHmHr'”‘**% ?  <« 


u#  .  v  ,  nyr.  Wytoc  I 
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For  the  changes  in  curvature,  we  shall  use  the  quantities  pre¬ 
viously  used  b>  Donnell  (8>  for  perfect  cylinder,  vs. 

x,  -  «■.,  x,  -  ^7.  xw  -  ^jjr  [isi 

The  worn  oont  H  D  b>  tbe  uniform  external  pressure  acting  on 
an  sides  <a  the  cj  Under  is  given  b>  the  product  of  the  pressure 
and  the  change  in  volume  of  the  cylinder.  Thus  there  follows 
(16) 

V.(W  +  nl-~V(W.  -uO 
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We  skal’  eirr  puf>  this  complicated  exprea*.on for  IV’ to  tne  J  !!ow- 
ing 
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I*  .»  aiso  convenient  *„  consider  the  V,  V  and  Ur  displacements 
to  be  made  up  of  two  parts 

V  -  C  -  II  ■ 


IF  -  fT  +  ir 


+  (w  +  w*),s  t-  (w  +  «*)*  j  (c)  j 

where  9,,  9§,  and  are  the  membrane  "tresses  which  would 
occur  in  a  perfectly  circular  cylinder  and  which  were  assumed  to 
be  constant  in  deriving  the  foregoing  equations. 

For  the  case  of  uniform  external  hydrostatic  pressure  applied 
on  all  sides  of  a  perfectly  circular  cylinder,  and  neglecting  the  de¬ 
flection  of  the  shell  between  supports,  tbeie  results 

■'lM1 

Substituting  the  Relations  120)  into  Equation  [19c]  yields 
0V4*  7  V’hr,,,,  +  pR  (u»  +  is#)., 

+  £,<«  1-  J  —  r  “  0  ...  {21 J 

Now  define  a  stress  function  of  the  total  membrane  stresses 
F,  such  that 
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Using  Equation  [14]  and  the  stress-strain  relations  it  is  then 
easy  to  show  that 

W.*LI, 

|  R  "  +  R>  R< 

,  2ir.,m.,  »*«>.. \ 

+  ~  R>  W  ]  ■  |Z" 

As  we  are  interested  here  in  the  linear  problem,  and  as  we  have 
previously  neglected  the  deflection  of  the  shell  between  supports. 
EquMion  (23)  reduces  to 

V‘F  -  -  ~  a.,  |24| 

Equat/ons  [19a],  11961  (2l|,and  (24)  together  with  the  appro¬ 
priate  boundary  conditions  define  the  problem. 
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